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INTRODUCTION 


The  light  emitted  from  missile  trails  represents  an  exceedingly 
complex  phenomenon.  A  large  variety  of  physical  theories  have  been  of¬ 
fered  to  explain  this  luminosity.  Since  particulate  matter  is  generally 
found  to  some  degree  in  missile  exhausts  an  understanding  of  the  inter¬ 
action  of  particulate  matter  with  the  atmosphere  may  provide  a  key  to 
the  disentanglement  of  the  particulate  matter  contribution.  The  thesis 
here  is  that  if  a  particular  zone  of  luminsoity  is  observed  to  move  in 
a  certain  fashion  which  is  inappropriate  to  the  movement  of  particulate 
matter  then  the  invoking  of  the  particulate  matter  mechanism  cannot  be 
accepted. 

In  an  effort  to  provide  a  means  for  checking  the  particulate  matter 
hypothesis  this  study  has  been  made.  It  attempts  to  systematically  obtain 
the  mathematical  formulae  for  the  velocity  and  distance  behavior  of 
particles  of  different  mass  and  initial  velocity  ejected  at  various 
heights  in  an  exponential  atmosphere.  Its  purpose  is  essentially  that 
of  an  engineering  type  study  for  application  to  the  above  cited  missile 
phenomena.  Hence  for  clarification  specific  examples  are  calculated  in 
detail  for  the  various  cases.  Also  a  more  extensive  tabulation  for 
horizontal  motion  With  drag  is  given  in  an  appendix. 

By  use  of  the  work  contained  herein  it  is  hoped  that  the  particulate 
matter  component  in  missile  trail  luminsoity  may  be  more  carefully  check¬ 
ed,  This  study  represents  a  small  effort  in  a  complex  over-all  program 
of  evaluating  the  missile  trail  radiation. 
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If  particulate  matter  is  released  in  the  upper  atmosphere,  at  a 
point,  with  a  common  initial  speed  in  all  directions,  the  drag  resistance 
of  the  atmosphere  will  slow  down  the  particles  in  such  a  manner  that  the 
space  loci  of  the  advancing  particles,  with  a  common  varying  speed,  will 
describe  successively  expanding  aspherical  surfaces.  That  is  to  say, 
the  surfaces  of  constant  speed  will  be  strongly  dependent  upon  the  vary¬ 
ing  density,  p(h),  of  the  upper  atmosphere,  which  is  essentially  expon¬ 
ential, 

p  =  p(h)  =  pQ  exp(-h/HQ)  (1) 

for  different  scale  heights  Hq. 

The  aerodynamic  drag,  D,  on  a  particle  of  speed  v  and  resisting 
area  A,  moving  through  such  a  density  field,  is  defined,  in  terms  of  a 
drag  coefficient  by  the  standard  formula 

(2) 

acting  in  such  a  manner  as  to  directly  oppose  the  motion  of  the  particle. 
The  force  of  gravity  can  be  taken  to  be  constant,  W  =  mg,  in  a  downward 
direction. 

The  combined  action  of  aerodynamic  drag  and  gravity  on  the  particle 
is  to  impart  to  the  particle  an  acceleration  a  defined  by  the  Newtonian 
equation  of  motion 
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ma  =  -  D  —  -  W  e,  (3) 

v  h 

where  the  velocity  vector  v  is  defined  by  the  equations 

a  =  dv/dt  (4) 


and 


v  = 


dr 

dt 


e 

r 


(5) 


where  e^  is  a  unit  vector  in  the  upward  direction  and  e^  is  a  unit  vector 
in  the  radial  direction,  out  from  the  axis  of  symmetry,  the  vertical  axis 
through  the  initial  release  point  r  =  0,  h  =  hQ . 

Equation  (3)  is  equivalent  to  the  two  componential  equations  of 
motion 


D  dr 
v  dt 


(6a) 


and 


2h  D  dh 
dt2  =  “  v  dt  " 


(6b) 


where 


v2  =  (dr/dt)2  +  (dh/dt)2 


(7) 
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The  initial  conditions  of  motion,  when  t  =  0,  are  defined  by  the 
statements: 


r  =  0,  dr/dt  =  v  sin  cc  , 
*  o  To’ 

h  =  h  ,  dh/dt  =  v  cos  cp  , 
o  o  To 


(8) 


the  parameter cpo  selects  a  specified  set  of  particles  issuing  from  the 
apex  of  a  cone  of  generating  angle  cpQ,  measured  from  the  vertical  line 
through  r  =  0,  h  =  hQ. 

If  the  times,  t,  of  observation  are  very  short  so  that  the  quantity 
gt  is  small  compared  to  v^,  one  may  neglect  the  gravitational  action  in 
describing  the  motion  of  the  particles.  In  this  event,  with  no  gravity, 
there  are  three  principal  cases  of  motion  with  the  following  differential 
equations  and  initial  conditions: 


Case  I:  Horizontal  Motion  with  Drag 


dvH/dt  = 


(9a) 


vH  -  dsH/dt 


(9b) 


v„  =  v  ,  sT  =  0,  when  t  =  0 
H  o  H 


(9c) 


where 


c  =  \  (CDA/m)  pQ  exp(-ho/HQ) 


(10) 
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Case  II:  Upward  Motion  with  Drag 


d^/dt  =  -  C  exp(-su/HQ)  v^ 


VU  =  d8U/dt>  8U  =  h  "  ho 


vIT  =  v  ,  sT  =  0,  when  t  =  0 
U  o  U 


(11a) 

(lib) 

(11c) 


Case  III:  Downward  Motion  with  Drag 


dvD/dt  =  -  C  exp(sD/HQ)  v^ 


vD  =  dsD/dt,  sD  =  hQ  -  h 


VD  =  Vo*  SD  =  °’  When  t  =  ° 


(12a) 

(12b) 

(12c) 


When  a  constant  gravity  field  is  introduced  into  the  description  of 
motion,  as  it  should  be,  two  additional  cases  may  be  added  to  our  sched¬ 
ule: 

Case  IV:  Upward  Motion  with  Drag  and  Gravity 


dVy/dt  -  -  C  exp(-sv/HQ)  vy  -  g 


vy  =  dsv/dt,  sv  =  h  -  ho 


Vy  —  0 ,  Sy  =  0,  when  t  =  0 


(13a) 

(13b) 

(13c) 
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v^?=  vo,  0,  when  t  =  0  (14c) 

An  analysis  of  these  five  cases  is  made  in  this  report.  The  sixth 
case  of  horizontal,  motion  (initial  release,  that  is)  with  drag  and 
gravity  is  not  discussed,  nor  is  the  general  problem,  described  by 
Equation  (3),  discussed.  In  order  to  illustrate  the  use  of  the  formulae, 
as  they  are  developed  in  each  successive  case,  application  is  made  to  a 
release  of  particles  at  71  Km,  with  an  initial  speed  vq  =  2500  cm/sec. 

The  numerical  examples  serve  as  a  means  of  comparing  the  various  features 
and  salient  points  of  each  case. 

In  the  formation  and  interpretation  of  missile  trails,  it  is  perti~ 
nent  to  know  at  what  altitude  one  can  expect  a  specified  drop  of  initial. 

velocity  in  a  specified  time,  and  to  know  the  penetrating  distances  s  , 

H. 

Sy,  and  Sy.  In  the  Appendix,  a  table  of  results  is  given  for  a  given 
set  of  initial  velocities  in  the  horizontal  direction,  for  Case  I, 

In  closing,  it  is  noted  that  this  paper  is  an  introductory  engineer" 
ing  study  only,  revealing  points  and  topics  worthy  of  further  study.  The 
fact  that  the  illustrations  used  are  couched  in  the  language  of  missile 
trails  should  not  obscure  the  evident  fact  of  the  generality  of  this 
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study  and  its  applicability  to  the  solution  of  problems  of  the  relative 
motion  of  bodies  in  an  exponential  atmosphere,  with  gravity. 
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HORIZONTAL  MOTION  WITH  DRAG 


In  this  section  we  consider  Case  I,  motion  in  a  horizontal  direction 
with  drag  only.  The  dynamic  equation  of  motion 


dvH/dt  =  -  C  vH  ,  vH  =  dsH/dt 


with  the  initial  condition 


v,  =  v  ,  when  t  =  0 
H  o’ 


by  direct  integration  yields  the  velocity  formula 


a/vK)  -  a/vo) .  c  tu 


which  can  be  written  in  the  forms 


v  /v„  =  1  +  C  v  t„ 
o  H  o  H 


(16a) 


c.-L-(i.i) 

Vh  V  } 


(16b) 


On  the  other  hand,  since  vu  =  ds„/d'  Equation  (9a)  can  be  written  as 

rl  ri 


dvR/dt  =  -  C  vH  (d8H/dt) 


which  has  the  integral 


vH  =  vq  exp(-C  sH) 
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since  v„  =  0,  when  s„  =  0,  That  is  to  say, 
ri  n 

vo/vH  =  exp(C  sH) 
so  that 

SH  =  C  in(vo/vH} 


(19) 


(20) 


By  using 

1  Vh 

C  -  (vo/vH)  -  1 


(21) 


from  Equation  (16),  we  can  express  su  in  the  significant  form 

ri 


H 


in(vo/vH) 
.(V0/VH}  "  1 


(22) 


Equation  (22)  defines  a  formula  for  the  horizontal  displacement,  s„,  of 
_  _ rl 


the  particle  in  the  time  t  =  tu,  when  the  velocity  ration  is  vu/v  . 

_ _ n  _  _ _ H.  O 


It  is  interesting  to  note  that,  for  a  fixed  value  of  VH/VQ>  this 
formula  for  states  that  the  displacement  s^  is  proportional  to  vQt, 
the  displacement  without  resistance,  when  C  =  0  (i.e.,  when  s  0,  say). 
Hence,  Equation  (22)  can  be  written  as 


SH 


H 


C=0 


in(vo/vH) 

(VoV  '  1 


(23) 
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this  ratio  is  illustrated  in  Figure  1.  Note  that  when  vu/v  =  1/2,  the 

displacement  with  drag  is  69%  of  the  displacement  without  drag,  i.e., 

8u / 8tt  =  697«;  when  vu/v  =  1/10,  then  8„/8u  =  25.57».  Lastly,  it 

H  Hc=o  H  °  H  Hc=o 

is  important  to  note  that  v„/v  can  have  any  value  from  unity  to  zero: 

ri  O 

from  Equation  (16a),  as  v  /v  -*•  0,  we  have  t„  -+  oo  ;  and,  from  Equation 

rl  O  rl 

(18),  as  v„/v  -+  0,  we  have  su  -*■  od  .  Practically,  of  course  this  is 

n  O  rl 

meaningless,  since  we  must  use  the  time,  t  ,  in  a  restricted  interval 

H. 

(since  we  neglect  the  effect  of  gravity). 

In  closure  of  this  discussion  on  horizontal  motion  with  drag,  notice 
that  we  can  find  the  altitude  hQ  required  to  have  sufficient  air  resist¬ 
ance  to  slow  the  particle  down  to  a  given  fraction  v.,/v  of  its  initial 
speed  in  a  given  specified  time  t^;  this  is  done  by  combining  the  two 
equations  for  C,  Equations  (1)  and  (16b),  namely, 

C  .  I  <CdA/„)  Po  exp(-ho/Ho)  =  ^  ^  ■  l)  «*> 

and  solving  for  the  density  function, 

p(hQ)  =  PQ  exP("h</H0)  (25a) 

in  terms  of  v  ,t  ,  and  v  /v„  in  the  form 

On  On 

P(V  -  tV  (  t;  - l)  <"/cda)  •  <25b> 

o  H  ^  H  / 
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HORIZONTAL  DISPLACEMENT  WITH  DRAG  SH 

HORIZONTAL  DISPLACEMENT  WITHOUT  DRAG  ’  SHC=0 


LTS2-3/0 


'I  to  too  v0/\ b 

1 1  I  I  I I I _ I _ I _ I  I  1  I  1  ■  J — I — I - 1 - 1 

/  ,5  o.t  'os  0.01  \fa/v0 

FRACTION  OF  INITIAL  HORIZONTAL  VELOCITY,  v0/vH 


Figure  1.  Variation  of  Horizontal  Displacement  with  Horizontal 
Velocity,  with  Drag  only  on  Particle. 
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The  right  hand  side  of  Equation  (25b)  is  known  as  soon  as  t^  and  vq/v^ 

are  selected  for  a  particular  release  vq  and  particle,  m/C^A;  hence 

p(hQ)  is  known.  To  find  hQ  one  can  use  Equation  (25a)  for  specified 

values  of  pQ  and  Hq,  which  can  be  found  in  a  standard  model  atmosphere 

handbook,  One  can  use  the  ARDC  Model  Atmosphere  tables  directly, 

entering  with  p  =  p(hQ)  and  reading  out  hQ,  the  altitude  corresponding. 

In  the  Appendix  will  be  found  results  of  calculations  for  particles  of 

radii  2| i,  ,5|i  and  .  03^i,  for  =  1  and  CD  =  2,  for  a  range  of  initial 

2  5 

velocities  v  ,  from  v  =  10  cm/ sec  to  v  =  5  x  10  cm/sec. 
o  o  o 

Equation  (25b)  defines  a  formula  for  the  atmospheric  density,  p(hQ), 
required  to  slow  a  particle  down  from  an  initial  speed  vq  to  the  value 
v..,  in  the  time  t„.  Equation  (25a)  can  be  used  to  determine  the  corre- 
sponding  altitude  hQ  (when  the  ARDC  Atmosphere  Model  tabulation  is  not 
available) . 
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II.  UPWARD  MOTION  WITH  DRAG 


In  this  section  we  consider  Case  II,  motion  in  an  upward  direction 
with  drag  only.  From  Equation  (11),  we  assert  that  Vy,  the  velocity 
upwards,  is  defined  by  the  differential  equation 

dVy/dt  =  -  C  exp(-Sy/HQ)  Vy2,  Vy  =  dSy/dt,  Sy  =  h  -  ho  (26a) 

with  the  initial  conditions 

v„  =  v  ,  8.  =  0,  when  t  =  0  .  (26b) 

U  o  U 

Since 

dWdt  =  (dVy/dSy)  (dSy/dt)  =  Vy  (dVy/dSy)  (27) 

we  can  express  Equation  (26a)  in  the  alternate  form, 

(l/Vy)  (dVy/dSy)  =  "  C  C  XP  (-Sy^)  ^  (28) 

which  has  the  integral 

£n(v  /v  )  =  -  C  H  [1  -  exp (-8  /H  )]  ,  (29) 

u  o  o  u  o 

Solving  for  the  displacement  Sy/H^,  in  Equation  (29),  we  see  that 


Since  the  left  hand  side  of  Equation  (30)  is  never  negative,  it  is  clear 

that  v  /vIT  is  restricted  to  the  range 

o  u 

1  -  <vft/vn)  £  exp  (OH  )  (31) 

o  u  o 

for  all  distances  Sy  ^  0.  In  other  words,  the  minimum  value  that  Vy/vQ 
can  attain  is  exp(-CHQ),  that  is, 

Wmin  =  exP(-CH0)  <32) 

When  Vy/vo  -*•  exp(-CHQ)  we  see  from  Equation  (30)  that  exp(-Sy/HQ)  ->  0, 

which  means  that  Sy/H^  -*  oo  .  As  the  particle  moves  upwards  it  is  slowed 
down  to  the  asymptotic  vy  =  vq  exp(-CHo)  as  it  recedes  away  from  the 
initial  point. 

Equation  (30)  defines  a  formula  for  the  vertical  accent  of  a  particle, 

8„,  in  terms  of  the  speed  ratio  v  /vIT  ^  exp(CH  ). 

U _  o  U _  o 

It  is  an  easy  matter  to  compare  the  displacements  in  the  horizontal, 

s„/H  ,  and  vertical,  s.,/H  ,  directions  under  the  condition  that  the  speeds 
Ho’  ’  U  o’ 

are  equal, 

VH/vo  =  VU/vo  “  exP<CH0)  •  <33> 

To  this  end  we  refer  to  Equation  (19)  for  vu/v  , 

ri  O 

vH/v  =  exp(=CsH)  =  vu/v0»  from  Equation  (33) 
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and  substitute  v  /vt1  =  exp(Cs„)  in  Equation  (30);  the  resulting  equation 

O  U  n 

is 

exp(-su/H0)  =  1  -  sh/Ho,  0  =  ®H  =  K0  *  <34> 

Equation  (34)  is  a  formula  relating  8„/H  and  S../H  for  the  condition 

HO  U  O 

that  vH/vo  =  Vu/vo, 

To  find  the  time  required  for  the  particle  to  slow  down  from  vq 
to  Vy,  in  ascending  the  distance  Sy  above  hQ,  we  turn  to  Equation  (29) 
and  write  it  in  the  form 

dsy/dt  =  vq  exp  f-[l  -  exp(-Sy/Ho)]  CH^  (35) 

which  can  be  solved  for  the  differential  of  time,  dt,  in  the  form 

dt  =  exp  (\l  -  exp(-Sy/Ho)]  CH^)  dSy  (36) 

Integration  yields  the  formula 

ty  =  ~  exp(CHQ)  F  U  exp[-CHo  exp(-z/H  )]  dz  (37) 

Equation  (37)  is  a  formula  for  the  time  t  =  ty  required  for  the  particle 
to  ascend  a  distance  Sy. 
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let 


To  find  ty  in  terms  of  vy/v0>  we  again  turn  to  Equation  (24)  and 


i  =  [1  >  exp (-8  /H )]  CH  =  Jn(v  /v„) 

U  O  O  O  u 


(38) 


so  that 


ds. 


ds 


d£/dt  =  C  exp(-su/Ho)  —  =  C[1 


dt 


U 


hence 


dSy/dt 


1  d  t 

C[1  -  (£/ CH  )]  dt  =  Vo  exP("^»  from  Equation  (35)  /  (39) 


We  rewrite  Equation  (39)  as 


Ho  exp  ( O 

dt  -  T  CHTt  C 

o  o  5 


(40) 


where  £  =  /n(v  /v  )  so  that  £  =  0  when  t  =  0  and  vTI  =  v  . 

o  u  u  o 


By  integration  of  Equation  (40)  we  obtain  the  formula 


H  rin(y  /v7T)  .v, 

=  —  T  °  cHXp.^  =  CHo  '  (41) 

o  'J  o  o 


U  v 


Equation  (41)  is  a  formula  for  the  time  t  =  t^  required  for  the  particle 


to  slow  down  from  v  to  vIT,  with  vTT/v  =  exp(  =  CH  )„ 
_  o _  U  U  o  r  o 
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We  now  take  up  a  brief  study  of  this  formula,  Equation  (41),  for 
ty  For  one  thing,  we  can  show  that  ty  -»  a>  as  vQ/Vy  exp(CHQ): 
from  Equation  (41)  (see  Figure  2) 


lim 


(W 


exp(CHQ) 


H 

t  =  lim  — 
U  n  V 

€  -»  0  O 


i  CH  -e 
o 


exp(x) 

CH  -x 
o 


dx 


=  lim 
€ 


H 


0  Vo 


.CH  -e  .  H 

o  dx  o  ,  , 

— -  =  —  lim 

CH  -x  v  n 

5  O  O  6  -+  0 


£n 


CH 


=  oo 


(42) 


That  is  to  say,  it  takes  an  infinite  time  for  the  particle  to  reach  its 
minimum  velocity,  vq  exp(“CHQ)  as  would  be  expected  intuitively. 

In  order  to  evaluate  t„  for  v„/v  =  exp(-CH  )  we  set 

U  U  o  r  o 

CH  -  x  =  y  (43) 


in  Equation  (44)  and  make  the  appropriate  reductions;  we  find  that 


t 


U 


H 

—  exp(CH  ) 
v  r  o 


i  CH 


CH  -B 
o  r 


exp(-y)  dy 
y 


(44) 


where 


P  =  £n(vQ/Vy),  Vy  =  vq  exp(-p),  0  1  p  <  CH, 


(45) 


If  we  now  let 


P  =  (1”4)  CHq,  0  =  (a=l 


(46) 
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so  that 


Jahnke  and  Ende  give  only  a  small  graph  of  the  function  -Ei(-x) ;  reading 
from  this  graph  one  can  construct  the  following  table  of  values  (Table  1) 


TABLE  1 


VALUES  OF  -Ei(-x) 


pCO 

-Ei(-x)  =  /  exp(-y)  y 

d  o 

dy 

X 

-Ei(-x)  x 

-Ei(-x) 

X 

-Ei(-x) 

0 

oo  ,  6 

.45 

1.2 

.15 

.1 

1.70  .7 

.40 

1.3 

.12 

.2 

1.20  .8 

.27 

1.4 

.11 

,3 

.90  .9 

.24 

1.5 

.10 

.4 

.70  1.0 

.20 

1.6 

.09 

.5 

.55  1.1 

.16 

the 

other  hand,  the  formula 

for  ty  in 

Equation 

(41)  can 

proximated  analytically  when 


£n(v  /v„)  «  CH 
o  U  o 


(51) 


for  then 


1 


r'^n(vo/vu)  x 

cu  =  c T  J  exp(x)  (1  +  ciT  + 

o  °  o  o 


(52) 


so  that 


1  f  \x=in(v  /v  ) 

ty  =  ^r(exp(x)  +  [exp(x)/CHQ]  (x-l)J  , 

o''  /  x=0 


Simplification  of  Equation  (53)  results  in  the  equation 


CU  “  Cv  [ <vo/vu)-l] [ 1-(1/CHq) ]  +  Cv  CH  £n(vo/vu) 

O  0  0 


(54) 


or 


CU  =  Cv"  t(vo/vU)_l]  +  cTT  CH~  ([(vo/vU)  +  (55) 

O  O  n  I 


Equation  (54)  is  a  formula  for  the  time  t  required  for  the  ascending 

particle  to  slow  down  from  v  to  v„,  v  /v„  «  exp(CH  ). 
r  o  U*  o  U  o 


If  we  now  compare  t  with  t  by  requiring  that  v  /v  =  v  /v  ,  it  is 
U  n  O  rl  O  U 


seen  that  Equation  (54)  is  equivalent  to 


t  _  t  +  — 
U  H  CH 


1-<vH/vo) 


[in(vQ/vH)-l  +  (VH/V0)1* 


Vo/vH  =  Vo/vU  =  exP<CV 


(56) 


It  is  clear  that  if  v  /v„  =  e  then  tTT  >  tu.  It  can  be  shown,  however, 

Orl  U  n 

that  t  >  t  even  when  1  =  v  /v„  =  e;  for  let  v  /vtI  =  exp[  l-|a(vu/v  )], 

Url  OH.  OH  HO 

0  <  4  <  1,  so  that  in(v  /v„)  =  l-n(v„/v  )  then  in  Equation  (56) 

OH  HO 


CU  tH  +  CH  . 

o 


vH/vo 


l-(vH/vo) 


(1-n)  >  t 


H 


(57) 
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In  other  words,  for  the  same  drop  in  speed  it  takes  a  longer  time  In 

ascending  than  in  moving  horizontally,  tIT  >  t„  when  v  /v„  =  v  /v„  $ 

U  H.  O  H  O  U  > 

exp(CHo).  This  is  of  course  to  be  expected  intuitively  and  serves  as 

a  check  on  our  solution. 


Equation  (54)  can  be  written  in  the  form 
ty  =  [l-(l/CHo)]  tH  +  (1/CHo)(vo/vh)(sh/Ho), 


(58) 


(v0/vH)  =  (vo/vU)  «  exp(CV 

in  order  to  clarify  the  previous  discussion  and  to  illustrate  the 
use  of  the  formulae,  we  consider 


Example  1 


We 


-5 


take  a  particle  of  radius  .5^=  5  x  10  cm,  with  mass 


-12  -9  2 

m  =  1.84  x  10  cm  and  resisting  area  A  =  7.85  x  10  cm  so  that 


(m/CpA)  =  1.2  x  10  ^  gm/cm^  with  Cp  =  2 
We  stipulate  that 


(59) 


<V0/VH°  -  <W  *  10' 


with  v  =  2.5  x  10  cm/sec 


(60) 


so  that  vu  =  vIT  =  250  cm/sec;  we  also  set  t  =  10  sec.  The  required 
n.  U  rl 

density  is 


22 


2 


x  9  x  (1.2  x  10'4) 


P(h0) 


V 

t 

_£  .  i 

ra 

v„ 

C^A 

H 

D  1 

2.5  x  104 


-8  3 

=  8.64  x  10  gm/cm 


(61) 


£ 

The  ARDC  Model  Atmosphere  tables  give  hQ  =  71  km  =  7.1  x  10  cm  as  the 
corresponding  attitude.  The  constant  C  has  the  value 


-8 

C  =  t-  (C  A/m)  p (h  )  =  --  —  ~  x  '10  =  3.6  x  10  4/cm 

20  °  2  1.2  x  10-4 


(62) 


To  find  the  scale  height  Hq  at  hQ  =  71  km  we  look  up  in  the  ARDC  manual 

-8  3 

and  select  two  adjacent  entries:  h^  =  71  km,  =  8.64  x  10  gm/cm  , 

-8  3 

h2  =  75.5  km,  p2  =  4.32  x  10  gm/cm  ;  we  then  use  p(h)  =  pQ  exp(-h/HQ) 
for  each  of  these  and  form  the  quotient 

(pl/p2)  =  exp[  (h2-h1)/HQ]  'f  (63) 

hence,  with  (p^/p2)  =  2,  (h2«h^)  =4.5  km,  we  have 

exp(4.5/Ho)  =  2  =  exp (.693);  Hq  =  6.49  km.  (64) 


With  H  =  6.49  km,  we  have 
o 

CHq  =  (3.6  x  10"4) (6.49  x  105)  =  233.64  .  (65) 

Since  in(v  /vtI)  =  £n(10)  =  2.303  «  233.64  =  CH  ,  we  have  (v  /v1T)  «  exp(CH  ) 
o  U  o  o  0  o 

so  that  it  is  legitimate  to  use  approximate  formulas.  We  need,  however, 

the  horizontal  distance  s  ;  we  use 

H. 
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(2.5  x  104)  =  6394  cm 


(66) 


s 


H 


M(vo/vh) 

2 . 302 \ 

L'w-1 

V  lt.  — 

o  H 

9 

from  Equation  (22).  Hence  in  Equation  (58) 


t 


U 


I  i 

v  H 

1  1 

r"  ch” 

\  o 

CH  +  ciT  vu  H  -  (1-00428>10  + 
o  H  o 

233. 64 j 

x  (10) 


6.394  x  103 
2.5  x  103 


which  works  out  to  be 

t  =  9.9572  +  .10947  =  10.06667  sec  (67) 


Further,  from  Equation  (34),  we  have 

exp(-s  /H  )  =  1  -  (s„/H  )  =  1  -  .009846  =  .99016  (68) 

U  O  ri  O 

so  that  sy/Ho  =  .009887  and  s^  =  6417.75  cm 

3 

In  summary:  if  v  =  2.5  x  10  cm/sec,  v  /v„  =  v  / vTT  =  10  in  t  =  10  sec, 

O  O  n  O  U  rl 

and  if  m/C^A  =  1.24  x  10  4  gm/cm2,  C  =  3.6  x  10  4/cm,  Hq  =  6.49  x  103  cm, 
then  p(h  )  =  8.64  x  10  3  gm/cm3,  h  =  7.1  x  108  cm,  su  =  6.394  x  103  cm, 

O  OH. 

3 

Sy  =  6.418  x  10  cm,  t^  =  10.0667  sec. 

Example  2  (Gegenbeispiel) 

We  now  take  up  an  illustration  in  which  our  formulas  fail:  with 
the  same  partice  as  in  Example  1  we  now  let  vq  =  5  x  103  cm/sec  and  put 
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the  same  velocity  drop  ratio  requirement  as  before,  Vy/v  =  Vy/v  =  1/10 

-10  1 

in  t  =  10  sec.  We  find  that  p(h  )  =  4.32  x  10  gm/cm  and  h  =  99  km. 
n  O  o 

The  constant  C  is  now  1.8  x  10  ^/cm;  the  scale  height  Hq  is  now 

4.48  x  10^  cm,  and  CH  =  .806.  Since  in(v  /vIT)  =  ^n(10)  =  2.302,  we 

o  o  U 

see  that  (v  /vy) y  exp(CHQ)  in  violation  of  the  minimum  velocity  drop 
ratio  restriction.  That  is  to  say,  at  the  height  hQ  =  99  km  it  is  im¬ 
possible  for  the  particle,  beginning  with  an  upward  velocity 
vq  =  5  x  10^  cm/sec  to  slow  down  to  a  speed  Vy  =  (1/10)vq  =  5  x  10^  cm/sec: 
the  minimum  value  of  (Vy/vQ)  attainable  is  exp(-CHQ)  =  exp(-.806)  =  .4466. 

At  such  a  high  altitude  (h  =  99  km)  and  with  such  an  initial  speed 

o  _ 

✓ 

(v  =  5  x  10  cm/sec),  the  drag  resistance  can  only  slow  it  down  to  the 

value  vIT  =  .4466  v  . 

_ U _ o 

If  gravity  acts  on  the  particle  it  will  halt  it  in  time  (if  vq  is 
not  an  escape  velocity).  In  the  next  section  we  bring  the  gravity  into 
play,  in  combined  action  with  the  drag,  for  particles  released  in  an 
upward  direction. 
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III.  UPWARD  MOTION  WITH  DRAG  AND  GRAVITY 


In  this  section  we  consider  Case  IV,  motion  of  a  particle  with  an 
initial  upward  velocity,  with  drag  and  gravity  combined,  acting  to  slow 
it  down.  The  upward  velocity  is  now  designated  by 

vy  =  (dsy/dt),  sy  =  h  -  hQ  ,  (69) 

The  dynamic  equation  of  motion  is  expressed  in  Equation  (13)  as 
(dVy/dt)  =  -  C  exp(-sv/HQ)  vy2  -  g  (70a) 

with  the  initial  conditions 

vy  =  vQ,  Sy  =  0,  when  t  =  0  .  (70b) 

To  integrate  Equation  (70a),  we  begin  by  replacing  dvy/dt  by  VydVy/dSy, 
and  multiplying  through  by  dSyj  we  find  that 

vV  dvV  =  "  C  vV2  exP(-8v/H0)  d8V  "  8  dsV  (71) 

or 

d(g  Sy  +  J  Vy2)  =  CHo  Vy2  d  [  tt*  (  »  Sy /Hq  )  ] 

which  can  be  written  as 

d(g  Sy  +  Vy2)  =  2CHQ(g  Sy  +  |  )  d  [  CX?  (  »  Sy  /Hq  )  ] 

-  2CHq  g  sv  d[exp(-.sv/Ho)]  (72) 
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In  order  to  simplify  the  appearance  of  Equation  (72)  we  let 


1  2 

u  =  8  8V  +  2  VV  (73a) 

w  =  exp(-sv/HQ),  dw  =  -  (w/Hq)  d  sv  (73b) 

we  now  have  Equation  (72)  in  the  form 

du  =  2CH  u  dw  -  2CH  g  s„  dw  (74) 

o  o  °  V 

The  relationship  dw  =  -  w  d( s^/HQj reduces  this  to 
du  =  -  2C  u  w  d  8^  +  2C  g  w  s^  ds^ 

which  can  be  written  in  the  standard  form  of  a  linear  first  order  ordinary 
differential  equation: 

(du/ dSy)  +  2C  w  u  =  2C  g  w  s  (75) 


An  integrating  factor  of  Equation  (75)  is 


u  =  exp 


2C  w  ds,.  =  exp(-2CH  w) 
V  r  o 


so  that 


=  2C  g  w  Sy  exp(-?CHQw) 


(76) 


(77) 


is  an  equivalent  canonical  form  of  Equation  (75).  Note  that  w  =  1  and 
2 

u  =  (l/2)v  when  s,r  =0,  (t  =  0),  so  that  by  integration  of  Equation  (77) 
o  v 
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we  obtain  an  integral  of  motion,  namely. 


u  exp(-2CH  w) 
o 


7  v  2  exp (-2CH) 
z  o  o 


=  2C  g 


exp (-2CHqw) 


d  sT 


(78) 


Note  that 


exp (-2CHqw) 


d  s. 


p  w 

H  /  s„  exp(-2CH  w)  dw 

o  J1  V  o 


2C  SV 


H 

exp(-2CHQw)  +  ~ 


exp(-2CH  x)  — 
o  x 


(79) 


so  that  Equation  (78)  reduces  to 


u  exp(-2CHQw)  =  - 


exp  (- 


2CH  )  + 
o' 


sv  exp(-2CHQw) 


+  g  H 


o 


exp(-2CHQx) 


dx 

x 


(80) 


or,  finally, 


u  =  —  vq  exp 


-2CH 


,d-w)] 


+  g  s, 


+  g  Hq  exp(2CHQw) 


exp (-2CHqx) 


dx 

x 


(81) 
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But  U  =  g  Sy 


+ 


1  2 
2  VV  ' 


so  that  Equation  (81)  implies  that 


exp  [-2CH  (l-wv)j  -  w  g  Hq  exp(2CHQwv)  f 


dv 

exp(-2CHox)~  (82a) 


where 

w 

w  =  wy  =  exp(-sv/HQ),  (dwv/dt)  =  -  jp  vy  (82b) 

Equation  (82)  defines  a  formula  for  the  velocity  of  ascent  under  the  com- 
bined  action  of  drag  and  gravity  in  terms  of  the  ascent  distance  s„. 


As  a  check  on  this  formula,  it  can  be  easily  shown  that  it  does 
satisfy  the  dynamic  equation  of  motion  Equation  (70a). 


With  gravity  and  drag  acting  on  the  particle,  in  its  ascent,  it  is 
possible  that  the  particle  will  slow  down  to  a  halt,  at  which  point 
Vy  =  0,  momentarily,  before  the  particle  then  begins  to  descend,  under 
the  action  of  gravity.  We  designate  this  vv  =  0  condition  by  a  sub 

asterisk  on  w^,  that  is  w  =  w^  when  v^  =  0.  From  Equation  (82a),  with 

« 

Vy  =  0,  we  see  that  w*  is  defined  by  the  integral  equation 


exp (-2CH  x)  —  =  -  °  • 
r  ox  2g  H 


exp (-2CHq), 


w*  =  exp(-s*/Ho) 


(83) 


2 

Our  problem  to  find  w^  in  terms  of  vq  /2gHQ  and  2  CHq.  In  Equation  (83) 
if  we  set  x  =  1-y,  our  definition  of  w^  becomes 
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" 1_w*  dv  vo 

exp(2CHoy)  fL  =  — 

o  o 


Note  that 


1  -  w*  =  1  -  exp(-S*/Ho)  =  s*/Ho  +  if  (s*/Ho)  «  1;  (85) 


so  that  y  is  defined  between  D  and  (s*/H  )  «  1,  i.e.,  if  condition 

”  o 

Equation  (85)  holds,  a  binomial  expression  is  valid  on  (1-y)  ^  in 
Equation  (84).  Indeed,  with  (s*/H  )  «  1,  Equation  (84)  becomes 

"  O 


.  s  . /H  +  . 
*  o 


O  ^  o 


exp(2CHQy)  (1  +  y  +  ...)  dy 


1  1  y=s*/H 

=  exp ( 2CH  y )  +  - r  exp(2CH  y)  (2CH  y-1)  °  (86) 

2  o  0  (2CH  )2  °  0  o 

o 


a  =  (1/2CH  ) ,  |.  =  2Cs*,  ag.  =  (s./H  )  «  1, 


Equation  (86)  can  be  written  as 


»o  2  2  2 

—  =  a  exp(|^)-l  +  a  g*  exp (g*)  -  a  exp(^)  +  a  , 


«£*  =  «  1 


The  inequality  ag^  «  1  permits  a  reduction  of  Equation  (88):  multiply 


Equation  (88)  through  by  g.  and  rearrange  terms,  in  the  manner 
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2 

V 

5*  2g ir  =  [a|*  +  (a|*)2]  exP(U>  -  al*  -  0/2  I*  exp(^)  +  a2  I*  (89) 
2 

and  since  (a|.v)  «  a|^  we  infer  that  it  is  valid  to  write  Equation  (89) 

as 

2 

v 

o  2  2 

=  a  exP^*)  -  a  "  a  exP  (!*)  +  a 

°  o 

=  a jexp (!*)-! j  -  a2jexp(|*)-l] 


or  as 


v 

2gH  -  [®xP(l*)”lj  > 


“I*  «  1 


(90) 


in  Equation  (90)  we  obtain  a  formula 

(91) 

and  since  =  208^,  a  =  (l/2)CHo>  we  infer  that 

(2CH  )2  v  2 

exp ^Cs^)  =  1  +  2CH  -1~  2gH  *  <<;  *  '  ^) 

o  6  o 

Equation  (92)  defines  a  formula  for  the  maximum  height,  =  s  ,  attained 

by  the  particle  in  its  ascent,  under  the  combined  action  of  gravity  and 

drag,  provided  (s./H  )  «  1. 

_  _  «  o  _ 


By  simply  rearrangement  of  terms 
for  exp (£*),. 

2 


exp(|*)  =  1  + 


a(l-a)  2gH 
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We  pause  at  this  point  to  consider 


Example  3 

3 

We  take  the  data  of  Example  1:  v  =  1.5  x  10  cm/sec, 

H  =  6.49  x  10^  cm,  C  =  3.6  x  10  \cm,  CH  =  233.6;  then  with 
o  o 

g  =  980.616  cm/sec2, 

v  2  (2CH  )2  v  2 

2gH  =  •00491>  2CH  -1  2gH  =  2.29924  ,  (93) 

o  o  °  o 

hence  in  Equation  (92) 

exP(2Cs*)  =  1  +  2.29924  =  3.2992  exp (1.194) 
so  that 

(s./H  )  =  (1 . 194/2CH  )  =  (1.194/467.2)  =  .002555  «  1  .  (94) 

w  o  o 

The  criterion  (s^/Ho)  «  1  is  amply  satisfied  so  that  it  is  legitimate 
to  use  the  formula  for  s*  in  Equation  (94).  Hence,  the  maximum  ascent 
distance  which  the  particle  can  attain  is 

s*  =  .002555  Hq  =  1658.19  cm  (95) 

3 

when  the  initial  upward  velocity  is  vq  =  2.5  x  10  cm/sec  and  the 

altitude  is  hQ  =  71  km,  under  the  combined  action  of  drag  and  gravity. 

If  gravity  alone  were  to  act  on  the  particle,  it  would  be  stopped  in 

t  =  v  /g  =  2500/980.616  =  2,549  sec  and  it  would  ascend  a  distance 
8  o 
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8  =  (l/2)v  t  =  3186.25  cm.  Hence  the  effect  of  drag  combined  with  the 
g  o  g 

gravity,  is  to  cut  the  distance  from  s  =  3186.25  cm  to  =  1658.19  cm, 

g 

or  to  reduce  the  free  ascent  by  1528.06  cm.  Recall  from  Example  1  that 
the  particle  ascends  a  distance  s^  =  6418  cm  in  t^  =  10.0667  sec  while 
reducing  the  velocity  from  vq  =  2500  cm/sec  to  v^.  =  250  cm/sec,  under 
the  action  of  drag  alone. 

Example  4  (Gegenbeispiel) 

We  take  the  data  of  Example  2:  vq  =  5  x  103  cm/sec, 

C  =  1.8  x  10  6/cm,  h  =99  km,  H  =  4.48  km,  CH  =  .806;  then 

o  *  o  *  o 

v  2  (2CH  )2  v  2 

i§r  -  284-543>  icrrr  iir=  1208-163 

o  o  o 


hence 


exP(2Cs*)  =  1  +  1208.163  =  1209.163  =  exp(7.0977) 
so  that 

(s*/Hq)  =  (7.0977/1.612)  =  4.403  >  1  )  (96) 

this  value  of  (s./H  )  =  4.403  is  a  violation  of  the  restriction  that 
w  o 

<s*/V  «  1:  Equation  (92)  cannot  be  used  in  Example  4. 

In  working  Example  3  we  did  not  consider  the  time  required  to 
reach  the  maximum  ascent.  This  points  up  the  need  to  determine  the  time 
ty  required  for  the  particle  to  slow  down  from  vq  to  v^  while  ascending 
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a  distance  Sy  under  the  combined  action  of  aerodynamic  drag  and  gravity. 
To  find  ty  we  need  dty  which  we  can  define  by  use  of  Equation  (82): 

dtV  =  -  ^  dwV«  w  =  exP(‘8v/Ho>  •  <97> 


At  least  theoretically,  we  have  ty  defined  by  the  integral  of  Equation 
(97),  namely, 


wv'=  exp(sv/Ho) 


(98) 


where  Vy  is  defined  by  Equation  (82a). 

Equation  (98)  defines  a  formula  for  the  particle  to  ascend  a  distance 
sv  against  the  combined  actions  of  aerodynamic  drag  and  gravity. 


On  the  other  hand  since  Vy  =  dSy/dt,  we  can  write  dt  =  dSy/Vy  and 
then  express  ty  in  the  form 


In  either  case,  by  using  Equations  (98)  or  (99),  it  is  necessary  to  have 
v.y  as  a  function  of  Sy  or  Wy  =  exp(Sy/Ho).  We,  therefore,  turn  our 
attention  to  the  reduction  of  the  integral  in  Equation  (82a): 
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/ 
w  W, 


1  pl-w  J 

exp(-2CHox)  ^  .  exp(-2CHQ)  /  V  exp(2CHQy)  (100) 

V  ° 


Note  that 


1  -  WV  -  1  '  exp(-sv/Ho)  =  sv/Ho  +  ...  if  sv/Ho  «  1  / 


so  that,  in  Equation  (100) 


I  =  exp(-2CHo) 


r 


s  /H  +  ... 


[exp(2CHQy)](l+y)  dy 


=  exp  (-1/a)  a(l-a)[exp(gv)-l]' 


(101) 


by  analogy  with  the  development  for  w^,  where 


a 


2CH 


»  ly  *  ^8y  »  aly  =  ^Sy^o^  ^  ^  * 


(102) 


note  that 


2CHo  wv  "  a  "  h  +  *  *  * » 


!  -  Wy  =a|v  + 


(103) 


2CHo(l-wv)  =  |v, 


2CH  w  -  -  = 
o  v  a 


-  + 


When  these  abreviations  are  used  in  Equation  (82a)  our  formula  for  v 


becomes 
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Set 


vy2  =  vo2  exp(-|v)  -  2gHQ  exp(-|v)  a(l-a)  [exp  (|y) -1  ] ,  1  (104) 

o 


2  2gHo  1 

'  *  —2  “<1-a)  ’  “  ■  2CT 

V  o 


(105) 


then  Equation  (104)  can  be  written  as 

(W2  =  exP(-|y)  [1-r2texp(|v)-l]]  ,  (sv/Ho)  «  1  ,  (106) 

Taking  a  square  root  of  both  sides  of  Equation  (106)  we  see  that 


VV  =  <dsv/dt)  =  vo  exP<-CV  ^(1  +  T2)  *  T2  exp(2Csv)  (107) 


so  that  by  simple  algebra 


dc  =  -  exp(Cs„) 

rv  V 

1  o 


M 


i  2 

l-nr 


dsv  »  (Sv/H0)  <<:  1 


2  exp(2Csy) 


and  by  integration 


t„  =  — 
V  yv 


exp(Cx) 


o  °  o 


dx,  (sv/Ho)  «  1 


1+Y 


-  exp(2Cx) 


(108) 


If  we  now  set  u  =  exp(Cx),  —  du  =  exp(Cx)dx,  then  Equation  (108) 

v/ 


reduces  to 
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— t 


v  cv  r 

0 


exp(Csv) 


du 


\ 


(sv/Ho)  «  1 


(109) 


N--u2 

2 

r 


which  integrates  at  once  to  yield  the  formula 


V  Cv  r 
o' 


sin 


rexp(Csv) 


\ji-nr2 


-  sin  ^ 

1.  JL._. 

sli+x^l- 

,  (8v/Ho)  «  1  (110) 


where 


2  2®Ho 


r  =  — i  (i-«). 


a  = 


2CH 


It  should  be  noted  that  Equation  (91)  can  be  written  as 


exP(2Cs  )  =  i  +  (i/r2)  =  (r2+D/r2 


so  that 


exp(Cs  ) 


r?±l 


rexp(Cs  ) 


or 


\f 


1+y 


(111) 


the  implication  is  that,  if  tTT  =  t.  when  s„  =  s.  and  v„  =  0,  then  in 

y  rC  V  a  y 

Equation  (110) 


t*  Cv  x 
o' 


it  .  -1 

2  "  sin 


Wi+r 


■  •  <VV  «  l-  <112> 
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Equation  (110)  defines  a  formula  for  the  time,  t^,  it  takes  for  a 

particle  to  ascend  a  distance  s^  with  a  velocity  drop  from  vq  to  v^, 

provided  s^/H^  «  1>  with  drag  and  gravity  opposing  the  motion  of  the 

particle.  Equation  (112)  defines  a  formula  for  the  least  time  required 
for  the  particle  to  reach  a  position  where  it  comes  to  a  halt  momentarily, 

v  =  0,  at  a  distance  s  above  the  starting  point. 

_ V _ _ _  _ _  " 


Example  5 


We  take  the  data  of  Examples  1  and  3:  h  =71  km,  v  =  2.5  x  10 

o  o 

cm/sec,  C  =  3.6  x  10  ^/cm,  H  =  6.49  x  10^  cm,  2CH  =  467.2,  sir  =  s.  = 

o  o  ’V* 

1658.19  cm,  exp(2Cs.)  =  3.2992,  exp(CsJ.)  =  1.81638,  (2gH  /v  2)  =  203.6528, 
a  =  (l/2CHo)  =  .0021404,  1-a  =  .99786,  a(l-a)  =  .0021358,  r2  =  .4349616, 
r  =  .65952,  cot_1y  =  .98756,  sjl +y2  =  1.1979,  r/\)l+T2  =  .550563  = 
sin(.583044lrj);  Cvj-  =  .593568,  (l/Cvor)  =  1.684726;  (T/Jl+T2)  exp(Cs^)  = 
(.550562)(1. 81632)  =  1  =  sin(l. 57079  lr) ) ;  hence 


t 


V 


=  1.684726  (1.57079  -  .583044)  =  1.6641  sec 


cot  V  _  .98756 

Cvr  '  .59357 
o 


We  conclude  that  it  takes  t^  =  1.6641  sec  for  the  particle  to 
ascend  to  its  maximum  displacement  s^  =  1658.19  cm,  where  the  particle 
momentarily  comes  to  rest  v^  =  0,  under  the  combined  action  of  drag  and 

3 

gravity,  beginning  with  a  velocity  vq  =  2.5  x  10  cm/sec  upwards.  As  a 
point  of  comparison,  to  ascend  a  distance  s^  =  1658.2  cm  with  an  initial 
velocity  vQ  =  2500  cm/sec  and  be  acted  upon  by  gravity  alone,  it  takes 
.7836  sec. 
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From  Equation 


We  now  let  vQ/Vy  =  10  and  seek  the  distance  s^. 
(106)  by  simple  algebra  we  derive  the  formula 


exp(Csv)  = 


1+r^ 


(vV/vo)2  +  ^ 


(sv/R0)  «  i 


(113) 


Equation  (113)  defines  a  formula  for  the  distance  of  ascent,  s^,  of  a 

particle  under  the  combined  action  of  drag  and  gravity,  in  terms  of  the 
velocity  ration  Vy/vQ,  provided  sv  «  V 

Thus  with  1+y2  =  1.434016,  (Vy/vQ)2  +  y2  =  .01  +  .4349616  =  .444916, 
we  have  exp(Cs^)  =  \j 3.225139  =  1.7958  =  exp(. 58545)  so  that  (s^/Hq)  = 

(.58545/CHo)  =  (.58545/233.6)  =  .0025062  <  1,  s  =  1626.52  cm;  that  is 
to  say:  the  particle  ascends  a  distance  s v  =  1626.52  cm  against  drag 
and  gravity  while  dropping  its  speed  from  vq  =  2500  cm/sec  to 
Vy  =  250  cm/sec.  To  find  the  corresponding  time,  t^,  we  use  Equation 
(114,  wherein  Cs^  =  .58545,  exp(Csv)  =  1.7958,  yexp(Csv) /J 1+y2  = 

.988701  =  sin(81  .379)  =  sin(1.4203  r  ),  (y/^l+y2)  =  .550563  = 
sin(.58304U)  ),  hence  t  =  1.684726  (1.4203  =  .58304)  =  1.4105  sec; 
it  takes  t  =  1.4105  sec  for  the  particle  to  lose  90%  of  its  initial 
speed,  to  go  from  vq  =  2500  cm/sec  to  v^  =  250  cm/sec,  working  against 
drag  and  gravity,  in  going  a  distance  s^  =  1626.52  cm. 
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IV.  DOWNWARD  MOTION  WITH  DRAG 

In  this  section  we  consider  Case  III,  motion  downwards  with  drag 
only.  We  designate  the  velocity  vector  downwards  by  v^  and  define  it 
by  the  differential  Equation  (12a), 

dvD/dt  =  -  C  exp(sD/HQ)  vD2  ,  vD  =  dsD/dt  ,  sD  =  hQ-  h  (114a) 

and  have  it  satisfy  the  initial  conditions 

vQ  =  vq  ,  sD  =  0  ,  when  t  =  0.  (114b) 

•Following  the  techniques  used  previously,  we  find  that 

(l/vD)(dvD/dsD)  =  -  C  exp(sD/Ho)  (115) 

is  an  equivalent  form  of  Equation  (114a)  and  that 

in(vD/vQ)  =  -  [exp(sD/HQ)-l]  CHq  (116a) 

and 

(W  =  exp  |-[exp(sD/HQ)-l]  CHq |  ’  (116b) 

are  first  integrals. 

The  equations  in  Equation  (116)  defined  a  relationship  between  the  velocity 
ratio  vp/vQ  and  the  descent  distance  s^  for  a  particle  descending  into  an 

exponential  atmosphere  under  the  action  of  drag  alone. 
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By  simple  rearrangement  of  terms  we  transform  Equation  (116)  into 


exP(8D/H0)  =  1  +  CH"  in(vo/vD) 
o 


(117) 


It  is  of  basic  interest  to  compare  the  values  of  the  displacements  s^, 
Sy,  and  sD  at  the  instants  their  corresponding  speeds  are  equal,  that 
is  when  v  /v„  =  v  /v  =  v  /v  <  exp(CH  ). 

O  H  O  U  O  D  O 


From  Equation  (19),  in(v  /v„)  =  Cs  ,  so  that  in  Equation  (117) 

O  rl  rl 


exp(sD/Ho)  =  1  +  ir 
o 


(118) 


since  exp(-s  /H  )  =  1  -  (s  /H  ),  from  Equation  (34),  it  follows  that 
U  O  rl  O 


exP ( SD /Hq )  =  2  -  exp(-8y/HQ) , 


(W  =  (vo/vD)  =  (vo/vU)  =  eXp(CHo} 


(119) 


In  Figure  3  is  shown  a  comparison  of  the  displacements  Sy/HQ  and  sd/Hq 
with  s  /H  for  the  condition  of  equal  speeds. 

H  O 

To  find  a  time  dependency  relationship,  we  begin  with  Equation  (116b), 


vD  =  (dsD/dt)  =  vq  exp(-[exp(sD/Ho)-l]  CHq 


(116b) 


and  express  the  element  of  time,  dt,  as 


dt  =  expj[exp(sD/Ho)-l]  ChJ  dsD 


(120) 
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which  can  be  integrated  directly,  yielding  the  formula 
tD  =  [exp(-CHo)/vQ]  f  D  exp [ exp (£/H  )  CHq]  d£ 

J  J  o  u 

To  find  an  alternate  formula  for  t^  we  let 
i  =  [exp(sD/HQ)-l]  CHq  =  in(vQ/vD) 

in  Equation  (116b),  so  that 

(d$/dt)  =  C  exp(sD/Ho)  —  =C  [lr  (£/CHq)]  ~  > 

hence,  by  solving  for  ds^/dt  in  Equation  (123)  and  using  the 
relationship  in  Equation  (116b),  one  finds  that 

(<Vdt)  *  c  TTT7ciT  df  *  vo  e*P(-£> 

which  implies  that 


dt  _  5°  exEiLL  dt; 

C  v  CH  +  5  ^ 

o  o 


and  that 


\  pt  n(v„/vn) 

tD  “  v  J 

o  o  o 


o  D  exp.fr). 

CH  +  y  y 

n 


since  £  =  0  when  t  =  0,  from  Equation  (122), 


(121) 

(122) 

(123) 

velocity 

(124) 

(125) 

(126) 
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Equations  (121)  and  (126)  are  formulas  for  the  time,  t^,  of  descent, 
against  aerodynamic  drag,  in  terms  of  the  distance  s^,  in  the  former, 
and  in  terms  of  the  velocity  ration  drop  vq/v0>  ln  the  latter. 


In  the  event  that  VQ/vp  «  exp(CHQ)  the  kernel  of  the  integral  in 
Equation  (126)  can  be  expanded  in  a  series  and  integrated  termwiee,  thusly, 


CD  =  C  fr  "  t  -  cT  cr  ;r 7; +  1  -  -r  » -T  <<  exp(cV  .  (127) 
Do  o  o  D  D  D  '  D 


Equation  (127)  is  a  formula  for  the  time,  t^,  required  for  the  particle 

to  slow  down  from  a  downward  velocity  v^  to  the  value  v^,  provided 

v  /v  «  exp(CH  ).  To  appraise  this  formula  suppose  that  v  /v^  =  v  /vTT 
°  u  o  O  U  O  H 

then,  from  Equation  (15),  the  first  term  on  the  right  hand  side  of 

Equation  (127)  is 


t  =i  i-.-L 

H  C  VH  Vo 


so  that 


tH  ”  Cv  CH  (vn  in  vn  +  1  "  vn  ’  vn  "  v  <-<  exP(CH0)  •  (128) 

o  o  \  D  D  D  D  H 


An  alternate  form  of  Equation  (128)  is 


CD  "  CH 


t„  V  V„\  V  V 

-  in - 1  +  —  ]  ,  —  =  —  «  exp(CH  ) 

V  V  v  I  V  V  r\ 

CH  1 - —  D  °  '  D  VH  . 

°l  Vo 


(129) 
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Example  6 


We  take  the  data  of  Example  1:  hQ  =  71  km,  vq  =  2.5  x  10^  cm/sec, 

5  -4 

H  =  6.49  x  10  cm,  C  =  3.6  x  10  /cm,  CH  =  233.6;  we  choose 
o  o 

v  /v-.  =  v  /v„  =  1/10,  £n(v  /v_)  =  2.30259,  su  =  6.394  x  10"^  cm, 

O  D  O  n  O  u  H. 

s„/H  =  .009852;  hence 

it  O 

exp(s_/H  )  =  1  +  (s„/H  )  =  1.009852  (130) 

so  that  s^/H^  =  .0098035  and  s^  =  6362.47  cm  (which  compares  with 
8^  =  6418  cm  and  s^  =  1628.6  cm)  is  the  distance  the  particle  descends 
against  drag  alone,  in  reducing  the  downward  velocity  from  vq  =  2500  cm/sec 
to  vD  =  250  cm/sec.  To  find  the  time,  t^,  we  use  Equation  (127):  with 
v  /v_,  =  v  /vu  we  have  t  =  10  sec  (see  Example  1),  or  else  note  that 
Cv  =  .9,  v  /vn  -1  =  9;  further  v  /v^  in(v  AO  =  20.72331,  so  that 

O  O  U  O  D  O  V 

t  =  10  -  .055761  =  9.9442  sec  .  (131) 

The  time  required  for  the  particle  to  descend  the  distance  s^  =  6362.47 
cm  while  going  from  vq  =  2500  cm/sec  to  v^  =  250  cm/sec,  is  tp=9.9442  sec 
(compare  with  t^  =  10  sec,  t^  =  10.0667  sec  for  the  condition  vQ/vp  = 


V.  DOWNWARD  MOTION  WITH  DRAG  AND  GRAVITY 

In  this  section  we  consider  Case  V,  motion  downwards  with  drag  and 
gravity.  In  many  respects  much  of  the  previous  analysis  is  followed 
anew;  however,  with  gravity  now  acting  to  accelerate  the  particle  while 
the  atmospheric  drag  tends  to  slow  it  down,  the  resulting  motion  has 
different  features  so  that  the  analysis  is  not  without  interest. 

We  designate  the  downward  velocity  by  v^and  define  it  by  the  dynamic 
equation  in  Equation  (14): 

dv^/dt  =  -  C  exp ( Sy /Hp  )  v^>+  g  (132a) 

where 

V  Vdt  •  h° ' h  <132b) 

with  the  initial  conditions 


v-^  v  ,  s„  =  0  when  t  =  0 . 

Jr  o  ^ 


(132c) 


The  differential  equation  for  in  Equation  (132),  differs  from  that 
for  Vy  in  Equation  (70a)  by  having  -Sy/HQ  replaced  by  and  -  g  by  g. 

We  are  led  by  analogy  to  replace  -Hq  by  Hq  and  -g  by  g  ^in  the  formula 
Equation  (82a)  for  Vy ^to  write  down  by  inspection  a  formula  for  v^: 


v^  =  vq^  exp[2CHo(l-w^J]  +  2gHQ  exp(-2CHQ  v^_)  j  *^"exp  (2CHqx)  (133a) 
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where 


w^=  exP(^/Ho)  ,  (dw^dt)  =  (w^Hq)^v  (133b) 

w  =  1  when  s  =  0  ,  t  =  0  ,  (133c) 

We  verify  that  Equation  (133a)  defines  a  velocity  v^  that  satisfies  the 
equation  of  motion  by  substitution. 


Also  since  s.  r*  1  as  t  -*•  0  it  is  obvious  that  v.  -+•  v  in  Equation 
(133a),  so  that  the  initial  condition  on  v  is  satisfied. 


Equation  (133a)  defines  a  formula  for  downward  velocity,  „  of  a  particle, 

subjected  to  the  combined  action  of  drag  and  gravity,  in  terms  of  the 

descent  distance  s*v 

_ 


To  find  the  time  t .  of  descent  we  can  use  either 


H 

dt  =  — —  dw^,  from  Equation  (133b), 


(134a) 


or 


(134b) 


as  definitions  of  an  increment  of  time,  dt,  and  by  integration  obtain  the 
formulas 
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t  =  H 

&  oK 


'  exp  W  d“^ 
0 


(135a) 


(135b) 


The  formulas  in  Equation  (135)  define  the  time  tj  -for  the  particle  to 
descend  a  distance  s^-against  drag  and  with  gravity;  the  function  v^,.is 
defined  by  Equation  (133a). 


When  x  is  replaced  by  1+y  in  the  integral  in  Equation  (133a),  the 


formula  for  v.  becomes 


vj  =  vo2  exP[-2CHo(w^l)]  G(w^ 


(136a) 


o  w 

C(^i  =  1  +  <28h0/v02)  J  1  exp(2CHQy)  rj-^  ,  w^_  =  expOs^/H^  (136b) 


When  Hq  the  function  w^becomes 


w^=  exp(s^Ho)  =  1  +  (^JHo)  +  . 


(137) 


so  that 


'^1  «P(2™oy)  ^  . 


s  /H  + 

o 


exp(2CHQy)[l-y  +  ...]  dy  (138) 
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If 


a  =  <1/2CH0)  ,  l  -  2Cs^,  ai  =  (s^/Ho)  «  1 

and 

x2  =  (2gH  /v  2)  a(i  +  a) 

o  o 

then  it  can  be  shown  that 

(^/vq)2  =  exp +  x2  [l-exp(-j^J,]  = 

X2  +  (1  -  X2)  exp(-j=),  (s^l^)  <<:  1 


(139) 


(140) 


That  is  to  say 


>  2  •>  2 


C^vQ)  -  x  =  (1  -  X  )  exp(2Cs^)  ,  (s^/Ho)  «  1 


(141) 


Equation  (141)  defines  a  relationship  between  the  downward  velocity,  v. 


&L 


and  the  descent  distance,  a,  for  a  particle  under  the  combined  action 


of  drag  and  gravity,  s^_<;<  Hp. 


To  find  the  time  t*,when  s  .  «  H  ,  we  set 

oA-  o’ 


u  =  exp 


(C^J 


(du/dt)  =  C  v 


(142) 


so  that 
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-- 1 


dt -L 
C  uv«^  Cv 


da 


cv»  JxV  ♦  a-x2)  ^  ^5?  " »! 

When  dt  is  integrated  we  obtain  two  formulas: 


du 


.  (143) 


\er~  Cv  x 

w  o 


-1 

fXexp  (CsJ)' 

-Inh-1  X  \ 

Sinn 

x/l-x2  1 

Wi ~2I. 

,  X2  <  1,  ^<<  1,  (144a) 

o 


„  1  . -1 
t/=  - — 7  cosh 

^  CVoX 


/Xexp 


(C^> 


v/-2 


X  -1 


-  cosh 


■1  X 


W'2 


X  —  1/ 


,  X2  >  1,  ^<<  1  (144b) 

o 


Equation  (144)  defines  a  formula  for  the  time,  t^ ,  required  for  the  particle 

to  descend  a  distance  .  wherein  s*X<  Hq,  when  the  downward  motion  of  the 

particle  is  subjected  to  the  restraint  of  drag  and  the  acceleration  of 
gravity . 

Notice  that  Equation  (141)  can  be  expressed  as 


exp(CgP  = 


\J 


1  "X  ^  o 

—r—  —  ■  when  X  <  1  and  —  >  X,  — —  «  1  (145a) 

‘2  Vo  o 


<2^V  -  X 


or  as 


exp 


(C^>  = 


X2  -  1 


X2  -  (v  /vQ)2 


when  X2  >  1  and  X, 

v_ 


H 


«  1 


(145b) 
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If  Of  «  1,  then  X2  =  (2gHQ/vo2)a  =  (g/CvQ2) ,  2CHq  »  1; 


farther,  recall  that 


C  =  —  (C_A/m)  p  exp(h  /H  )  from  Equation  (10). 
z  u  o  o  o 


The  results  in  Equation  (145)  imply  that  when  «  1»  If  X  <1  then 

2 

the  minimum  value  of  v./v  is  X;  on  the  other  hand,  if  X  >1,  then 

0 

can  any  value  in  the  range  (0,1),  i.e.,  0  <  ^Vo  <  1‘ 


The  equation  in  Equation  (145)  defines  the  descent  distance,  in  terms 
of  the  velocity  ratloi  j^/vq>  provided  «  1. 


An  anomalous  behavior  of  the  velocity  vft / v  becomes  apparent  upon 

eC""  o 


examination  of  Equation  (141): 

v^/v6  =  X  \jl  +  (a-X^/X2)  exp(-2Cs^>  ,  (^/Hq)  « 


(147) 


if 


2  i  „  2 


( 1 1-X  |/X  )  exp(-2C^>1s  2f/X  ,  such  that  2f/X  «  1 


(148) 


then 


^Jvq  =  X  +  f  +  ....  when  exp(-2C^J)  <  (2fX/ 1 1-X2 1 ) ,  {^<<  1.  (149) 
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In  other  words 


v  =  4 g/C  ,  when  ^  ^ 


=  <  «  1 


(150) 


The  velocity  of  descent,  remains  essentially  at  the  constant  value 
4 g/C  for  a  finite  segment  of  its  descent,  from  an  altitude  hQ  with  an 
initial  velocity  vq;  the  parameter  C  varies  of  course  with  hQ  since 


c  =  j  (cDA/“0  P0  exp(-ho/Ho), 


In  this  connection,  if  we  assume  that  dj^dt  =  0  for  a  given  range  of 
values  s^,  we  can  set  d\g^/dt  in  Equation  (132),  the  dynamic  equation  of 
motion;  we  find  that 


v^,  =  (g/C)  exp(*^/Ho)  ,  (dv^/dt)  =  0 


(151) 


now  if  s i/H  «  1  then  exp (-s.  /H  )  ~  1  and  Equation  (151)  gives 
°  O 

£=  4&/c  in  agreement  with  Equation  (150).  Indeed  if  one  states  that 

dv../dt  =  -  5'  and  that  exp(s„/H  )  =  1  -  6"  ,  then 

&  jtr  ° 

-  5 '  =  -  C(l-S")  v^,  +  g  ,  so  that ,  if  v^=  4 g/C  +  e  ,  (152) 


then 


^7c 


=  \fg+5^/C(l-6")  -  Jg/C  -  nCT 


1+8 1 / g 
1-6" 


(153) 
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that  is , 


=  (1  +  \  B'/g)(l  -  \  6")  -1  =  7(7-'  6") 

Vg/c  Z  8  ’ 


(154) 


Equation  (154)  expresses  the  percent  error  in  assuming  the  velocity, 
a  constant ,  VgTc",  in  terms  of  the  variation,  in  the  rate  of 
change  of  v^  and  of  the  variation  6"  of  the  descent.  As  a  closure  to 
this  discussion  on  the  behavior  of  the  particle  for  constant  value  of 
v one  can  look  at  the  situation  from  a  physical  point  of  view:  the 
change  in  potential  energy  of  the  particle  goes  into  friction,  into 
heating  the  atmosphere,  by  the  frictional  work  done;  that  is  to  say, 
we  have 


W  ” c  exp(<^V  s4,lt. 

aJMnin 


(155) 


so  that 


2 


const 


CH 


^Tin 


exp(s^/Ho)  -  exP(<yHomln) 


c  >  exP^4H0)  «  1  (156) 


agreeing  with  the  previous  values  of  constant. 


Example  7 


We  take  the  data  of  Examples  1,  3,  and  5:  vq  =  2,5  x  10  cm/sec, 

h  =  71  km,  C  =  3.6  x  10'4/cm,  H  =  6,49  x  105  cm,  (1/a)  =  2CH  =  467.2, 
000 
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2gHQ/vo2  =  203.653,  a  =  .0021404,  X2  =  .43683,  X  =  .66093,  (XA/l-X2)  = 

.88072,  (XAT L-X?)  =  slnh  (.79447),  (1/C^>)  =  1.681127.  In  this  instance 
2 

X  <  1  so  that  we  must  choose  >  X  =  .6609;  but  this  is  not  as  re¬ 

vealing  as  assigning  a  value  to  t^J.n  Equation  (144a)  and  finding  the 
corresponding  value  of  To  this  end  we  take  t  ^  s  10  sec  and  solve 


for  s„ : 


sinh  *  |  — - —  exp(Cs)  =  Cv  X  t-  +  sinh”1 
•o*—'  o  otr 


\ji-x2 


\/ux2' 


=  6.74285  (157) 


so  that  exp 
hence 


=  480.24398  =  exp (13 .08205) : 


s<  =  (13.08205/3.6  x  10+4)  =  3.6339  x  104  cm,  (s  /H  )  =  . 05599. «  1 
cCr  Ob'  o 


We  infer  that  it  is  legitimate  to  use  these  approximation  formulas;  we 

4 

conclude  that  in  t  =  10  sec,  the  particle  will  descend  3.6339  x  10  cm, 
oCr 

assisted  by  gravity  but  opposed  by  aerodynamic  drag.  Compare  this  dis¬ 
placement  with  the  distance  the  particle  would  go  in  10  sec  under  3ravity 

4 

action  alone,  namely,  7.40308  x  10  cm.  The  speed  v  at  tA  =  10  sec  is 

rfr'  or' 

defined  by  Equation  (141)  wherein 


C^/v0)2  =  X2  +  (1-X2)  exp(-2Cs^)  = 


.4368314  +  [. 5631686/(480. 244)2]  =  .436831  +  ...  (159) 


(158) 
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so  that  0^«/vo)  =  X  =  ,66093,  v^=  1652.3  cm/sec.  Equation  (159)  in¬ 
dicates  that^yvQ  is  now  insensitive  to  the  value  of  8^ £<  Hq.  Notice 

that  if  we  take  s«^  =  6,748.9  cm,  then  exp(-Cs-^,  )=  exp(-2.4296)  = 

“^min  ^in 

(.8807196/10),  so  that  exp  (-2Cs-/  )  =  [( .880719)2/10o]  =  (1/100)  (X2/(l-X 

‘”min  L 

=  ^(101/100)  (X2/(l-X2J)j  -  (X2/(l-X2)) ,  hence,  by  reference  to  Equation  (141) 

where  we  have  exp  (-2Cs^)  =  ^(v^  /.v^)2/(l-X2j|  -  (X2/(l->^),  we  see  that  an 

error  of  less  than  10%  in  v„/y  occurs  when  s . ^  =  6748.9  cm  is  used 

°  ‘’min 

instead  of  4*^3=  36.339  cm  £  5.3  s^/  •  A  similar  discussion  reveals 
***  tnin 

4 

that  s n^,  =  (1/10)H  =  6.49  x  10  cm  causes  a  variation  in  the  value  of 

oCr  O 

max 

v  /vq  less  than  10%.  Hence,  for  ^jn  the  range  .0675  km  <  s^<, 

.649  km  the  variation  of  v-  =  1652.3  cm/sec  is  less  than  10%,.  Equation 

o<r' 

(144a)  is  used  to  determine  tw  ,  the  time  to  descend  to  s ^  ^  = 

°nvLn  min 

6748.9  cm  beginning  at  hQ  =  71  km: 


t^  =  1.0811272  [sinh_1(10)  -  .79447]  =  3.7048  sec 

min 


(160) 


since  (X/Jl-X2)  exp(Cs-^  )  =  10;  further  for  t*  ^  corresponding  to 

inin  'max 

s.  ^  =  6.49  x  10^  cm,  with  (X/xjl-X2)  exp(Csw  )  =  .88072  exp(23.364) 

^max  rnax 


exp (23, 237)  =  sinh(23.930) ;  hence 
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=  1.681127  (23.930  -  .79447)  =  38.894  sec  (161) 

^max 


The  time  period  of  constant  velocity  v*=  1652.3  cm/sec  persists  from 

cr 

=  3.705  sec  to  t^=  38.894  sec,  or  for  a  period  of  35.189  sec;  during 
35.189  sec  of  descent  at  a  constant  speed  1652.3  cm/sec  the  particle  falls 


a  distance  of  58,143  cm  (notice  that  when  we  compute  a 


^ 

max  min 


we 


obtain  64,900  -  6,749  =  58,151  cm.'.’). 
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